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Abstract 

This paper derives the anisotropic superexchange interactions from a Hub- 
bard model for excitations within the copper 3d band and the oxygen 2p 
band of the undoped insulating cuprates. We extend the recent calculation 
of Yildirim et al. [Phys. Rev. B VV, pp, 1995] in order to include the most 
general on-site Coulomb interactions (including those which involve more 
than two orbitals) when two holes occupy the same site. Our general results 
apply when the oxygen ions surrounding the copper ions form an octahe- 
dron which has tetragonal symmetry (but may be rotated as in lanthanum 
cuprate). For the tetragonal cuprates we obtain an easy-plane anisotropy in 
good agreement with experimental values. We predict the magnitude of the 

small in-plane anisotropy gap in the spin-wave spectrum of YB^CusOg. 
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I. INTRODUCTION 



The magnetic anisotropics of the family of compounds with structures similar to that 
of La2Cu04 (LCO) has been a subject of much current interest.llhi These materials are 
antiferromagnets,! with a dominant easy-plane anisotropy.0 The parent compound, LCO, 
is orthorhombic in the regime of interest. In this structure the oxygen octahedra surrounding 
each copper rotate by a small angle relative to their orientation in the tetragonal phase. This 
phase was found to have an antisymmetric Dzyaloshinskii-Moriya interaction, as allowed by 
the lack of inversion symmetry about the center of a Cu-Cu bond.0 More recently, members 
of this family which remain tetragonal, such as Sr2CuO2Cl2,0 have been studied. In many 
such compounds, for example, Sr 2 CuO2Cl2,0 Nd2CuO4,0 and Pr 2 Cu04lll the gap in the 
spin-wave spectrum at zero wave vector due to the easy-plane anisotropy has been found 
to be about 5meV, just as in LCO.0 

Since these systems involve copper ions in a 3d 9 configuration which have spin 1/2, the 
usually dominant mechanism of single-ion anisotropy does not come into play. Instead, 
anisotropy must be due to the anisotropy of the superexchange interaction. Microscopic 
derivations of the anisotropy energies have been carried out0~i on the basis of Anderson's 
theoryEl of kinetic superexchange, and Moriya's extension^ to incorporate spin-orbit inter- 
actions. These have been mostly confined to the consideration of the orthorhombic phase 
and were based on terms requiring the existence of a distortion. As a result, these calcula- 
tions produce anisotropy energies which are proportional to the distortion angle and which 
therefore vanish in the tetragonal phase. However, since the experiments cited aboveEUll 
indicate that the magnitude of the easy-plane anisotropy is independent of the struc- 
tural distortion, these theories, while giving correct information about the Dzyaloshinskii- 
Moriya interactionJlHl! did not provide a satisfactory basis for understanding the easy-plane 
anisotropy. Recognizing this fact, Yildirim et al.10 undertook an investigation of a model 
designed to calculate the anisotropy of the superexchange interaction so as to account cor- 
rectly for the tetragonal symmetry of the lattice. In their treatment all five crystal field 
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states of the copper 3d band (with a single hole on each copper ion) were taken into ac- 
count, as well as all the (occupied) 2p states of the oxygen ions. As a result, the exchange 
interaction associated with each Cu-Cu bond had biaxial exchange anisotropy. (I. e. all 
three diagonal components of the exchange tensor were in general different.) In that work, 
a model of Coulomb interactions was used which contained more terms than usual, but 
was still not completely general. Here we carry out the calculation with completely general 
on-site Coulomb interactions necessary to treat excited states with two holes on a single 
Cu ion, or on a single oxygen ion, assuming tetragonal site symmetry. As is known,0 such 
interactions can be parametrized in terms of only three parameters, the Racah parameters, 
A, B, and C, with A>B and A ^> C. In Refs. || and [7| it was shown that for tetragonal 
symmetry the exchange anisotropy vanishes if B = C = 0. Here we obtain results which 
include all contributions to the exchange anisotropy correct to first order in both B and C . 
Compared to the previous work, we find an additional contribution which slightly increases 
the biaxiality of the anisotropy, without sensibly changing the easy-plane anisotropy, which 
still agrees quite well with experimental results. Under certain approximations, our results 
can be applied to the orthorhombic phase of LCO. We should emphasize, however, that our 
aim in this paper is to calculate only the anisotropy of the exchange interaction. This point 
is discussed at the end of Sec. III. 

Briefly, this paper is organized as follows. In Sec. II we describe the model for which we 
calculate the exchange anisotropy. In Sec. Ill we describe the perturbative calculations, Sec. 
IV contains specific results for tetragonal symmetry, and in Sec. V we discuss and briefly 
summarize our results. We will attempt to make this paper self-contained, but readers 
wishing more details on the approach or the history of this problem are advised to consult 
Ref. |7| (YHAE). That reference also contains the details of the spin-wave spectrum which 
results from a calculation of the exchange anisotropy. 
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II. THE MODEL 



We start by considering the ground state of the Cu0 2 plane, when the kinetic energy 
is completely neglected. In that case, there is a single 3d hole on each copper ion and 
the oxygen 2p band is completely full. Since the spin of the Cu hole is arbitrary, this 
ground state is 2 7V -fold degenerate. The kinetic superexchange interaction is obtained as 
the effective interaction within this degenerate manifold when the kinetic energy is treated 
perturbatively. For this purpose we write the Hamiltonian as 



Here, 7icu C^o) describes the Hamiltonian of the copper (oxygen) ions and Hhop is the 
kinetic energy (hopping between Cu and O ions). To describe the Cu ions we work in 
a representation in which the crystal field Hamiltonian is diagonal. Diagonalization of the 
crystal field potential yields five spatial d-states, denoted by | a), with site energies e a . Each 
such energy is doubly degenerate due to spin a — ±1. In the systems under consideration, 
the five d-states | a) are determined by the tetragonal symmetry. Even in the orthorhombic 
phase the site symmetry of the copper ions can be taken to be tetragonal, since the crystal 
field is primarily generated by the neighboring oxygen ions, which, to a good approximation, 
still form an octahedron, albeit a rotated one. For tetragonal symmetry we label these crystal 
field states as | 0) ~ x 2 — y 2 , |1) ~ 3z 2 — r 2 , \ z) ~ xy, \ x) ~ yz, and | y) ~ zx, where 
the z-axis is perpendicular to the Cu02 plane and | 0) is the lowest energy single-particle 
state. Thus 



where d iaa creates a hole in the crystal field state a at site i with spin a. Here the first term 
is the crystal field Hamiltonian. The second term is the spin-orbit interaction, where A is 
the spin-orbit coupling constant and L Q/3 denotes the matrix elements of the orbital angular 
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T~t — T^cu + 7~(-o + T^hop • 



(1) 




(2) 



momentum vector between the crystal field states a and (3. The last term is the Coulomb 
interaction, where 

Uafrs = dri dr 2 ip a {ri)i! l3 (r2)—ip-y{r 2 )ip5(r 1 ) (3a) 

J J 7™i2 

= (a5\/3j) . (3b) 

in the notation^ of Table A26 of Ref. |2(| The second term in Eq. (|l|) gives the Hamiltonian 
of the oxygen ions. We assume that the spin-orbit interaction on the oxygen is much smaller 
than that on the copper and may be neglected. Hence we write 

= E e nPlnaPkna 
kna 

+ ^ E E (4) 

Z nin 2 n 3 n 4 aa 'k 

in which p\ na creates a hole in one of the three p-orbitals (denoted by n) on the oxygen at 
site k. Here the Coulomb matrix element is obtained analogously to Eqs. (Q). Finally, 7^hop 
describes the hopping between two neighboring oxygen and copper ions: 

Whop = E E tandiaaPkna + he , (5) 
iaa kna 

in which t^ n is the hopping matrix element and he denotes the Hermitian conjugate of the 
preceding terms. As mentioned, it is 7^ho P that lifts the degeneracy of the 2 7V -fold degenerate 
ground state. 

In order to derive the effective magnetic Hamiltonian it is convenient to start from the 
unperturbed Hamiltonian which contains the spin-orbit interaction exactly. To achieve this, 
we introduce the unitary transformation which diagonalizes the single-particle part of 7Ycu- 
As noted previously,!'@ the spin-dependence of this transformation is fixed by tetragonal 
symmetry, so that we can write 

[(*(<*)] (6) 

aa' 

where Cj acr / destroys a hole in the exact eigenstate of the Hamiltonian which consists of the 
crystal field and spin-orbit interactions. These states have a site label i, a state label (for 



which we use roman letters), and a pseudo-spin index a'. In Eq. (^) we define a(a) for each 
crystal field state | a) as follows: a(a) = a a is the Pauli matrix for a = x,y,z and a(a) = I 
is the unit matrix, for a = 0, 1. The 5x5 matrix m is the solution to 

£ a rn ab + taping = E b m ab (7) 

P 

with J2a m aa m oib = Sab, where 5 is the Kronecker delta function. Here C a p are related to 
the matrix elements of the orbital angular momentum vector and are given by 

Cq z = —iX, Cq x = Cq v = i\/2, 
Ci x = -Ci y = zV3A/2, 

C zx ^zy ^xy A/2. (8) 

Matrix elements not listed and not obtainable using C a/ 3 = C*^ a are zero. When A — ► 0, 
each state | a) approaches one of the states | a). Using this identification, the indices a also 
run over the values 0, 1, z, x, and y. We use the values of the parameters which are listed in 
Table I and are discussed in detail in YHAE. 

III. PERTURBATION EXPANSION 

We now divide the total Hamiltonian 7i into an unperturbed part, Tio, and a perturbation 
term Tii. The part 7io contains the single-particle Hamiltonians on the coppers and on the 
oxygens, and the leading on-site Coulomb potentials, 

7~(-0 — J]] E a C iacr Ci acr H — c iaa C iba' ' C iba' ' c iaa 

iau iab 
(J(T f 

+ 2 e nPknaPkna + PknaPkn'cr'Pkn'a'Pkna ■ (9) 

kncr hnn' 

For tetragonal site symmetry, we choose U = U aaaa = A + AB + C and Up = U nnnn . The 
perturbation Hamiltonian includes the kinetic energy and the remaining Coulomb interac- 
tions. Because of the transformation (^j), the hopping becomes spin-dependent. Thus we 
have 



n x = n hop + AH C , (io) 



in which 



^hop = EE ^o'Aac'Vkno + he, (11) 

with the 2x2 matrix 

^ = E^X«), (12) 

and A?ic describes the perturbation parts of the on-site Coulomb potentials, 



1 
2 

iabcd 



AH c = ^ AU ss' Sl s[(abcd)cl as cl bsf c icsi c ids/i 

+ 2 ^^ni™2ri3™4Plnio-Plra2o-'^ fera 30-'Pfcn40-5 (13) 

Jcctct' 



with 



AU SS , S1 seabed) = ^ AC/ Q/ 3 T(5 
a/375 

xm»#MK ft ) ff (*)L4 [^(/5)^(7)]-'«i- (I 4 ) 

Note that AU involves only the small Racah coefficients, B and C. 

The effective magnetic interaction 7i(i,j) between magnetic ions % and j is found by 
perturbation expansion with respect to Hi- All the perturbation contributions to 
should involve products of matrix elements which begin and end within the 2 iV -fold degen- 
erate ground-state manifold of Tlo, each state of which has one hole at each copper site, with 
arbitrary spin a. Denoting such a ground state by | ip ), and concentrating on perturbation 
terms which involve only two coppers i and j and one oxygen between them, it is convenient 
to replace | ipo) by c|ocr c ]oo-i c iOo-iQoo- I V'o)- This ensures that | ipo) indeed has exactly 
one hole on each copper i and j. Clearly the lowest-order contributions to the energy are 
of order t A . There are two possible channels in this order, which we denote by a and b. In 
channel a, the hole is transferred from one of the coppers to the oxygen, then to the second 



7 



copper. Afterwards, one of the holes on this second copper returns to the empty copper. 
Hence in this channel there are two holes on the copper in an intermediate state. In channel 
b, the hole is transferred from one of the coppers to the oxygen, and then a second hole is 
taken from another copper to the same oxygen. Then the two holes hop back, each to one 
of the initial coppers. Thus in channel b there are two holes on the oxygen in an intermedi- 
ate state. When the perturbation contributions coming from the Coulomb potential ATic 
[Eqs. (|10|), (|13D and flliD l are included, the on-site interactions on the copper are effective 
in channel a, and those on the oxygen appear in channel b. These contributions are of order 
t 4 AH c . 

The perturbation contributions to order t 4 are the same as those given by YHAE and 



arc 



HM(i,j) = £ g nn ,Tr{ff ■ Sdlt^a • S/£,^ } , (15) 

nn' 

in which 

S * = j 51 C i0a^)aa'Ci0a', (16) 



(7(7 



is the spin on the copper at site i in the orbital state | 0), t l Q n is the 2x2 matrix given by 
Eq. flTJ) and 

g nn ' = w + 7r — — - — (- + — ) • (17) 

The first term in g nn i results from channel a while the second arises from channel b. In 
deriving expression (|15D we have assumed that the single-particle energy E a=0 on the copper 
is equal to zero. 

From the form of T-C^(i,j) it is clear that it may yield anisotropic magnetic interactions 
only when the effective hopping between the coppers, t lk t k i , involves spin-flips. When the 
hopping conserves the spin, the t products are proportional to the 2x2 unit matrix and 
the trace in Eq. ( [15] ) gives a result proportional to Sj • Sj. This is indeed the case for 



tetragonal symmetry, as we discuss in the next sectionBQ In the orthorhombic phase of 
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La2Cu04, however, the effective hopping between the coppers is accompanied by spin-flip, 
and consequently T~C^ 1 \i,j) includes the antisymmetric Dzyaloshinskii-Moriya interaction, 
as well as symmetric magnetic anisotropics.^ 

We now turn to the contributions of the Coulomb potential. In particular, we consider 
contributions to the magnetic energy which are of order t A and first order in ATic- Following 
the approach of YHAE one finds that the processes from channel a (indicated by a subscript 
"a") yield 

Allans 



E 



,_, (U a + E a )(U + E b ) 

7ik Tkj 7jk 7ki 



/ lik 7k] ,jk ,ki 



m ab mp Q m l0 m Sa 



pk Zki 



x 

n' 6 ™' 



(18) 



where % <-> j denotes the sum of previous terms with i and j interchanged. Here it is 
convenient to classify the Coulomb matrix elements into two classes, by writing 

AUa^s = U^ lS + , (19) 

where U^l s is nonzero only if 

a(a) = a{5) , and a{(3) = a{^) , (20a) 

(2) 

while U^p^s i s nonzero only if 

a(a)cr(5) = Cicr^ = C 2 cr(/3)cr(j) , for fi = x,y,z, (20b) 

where C\ and C 2 are constants which may be imaginary. This classification will be of 
immediate use. For the elements denoted U^ s , the products a(a)a(5) and a{j3)a{^) are 
both proportional to the unit matrix. Therefore the first term in the square brackets of 
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Eq. ([Tj|) vanishes, and one is left with the second term alone. For the matrix elements 
denoted by U^ s , both products a(a)a(5) and a{(3)a{^i) are proportional to cx^, /i = x, y, z. 
In treating this contribution it is convenient to use the identity 

a^a ■ Sa^ = 2a fl S fM - a ■ S. (21) 

We see that the matrix elements s give rise to magnetic interactions which depend on 



the Cartesian index fi. Rearranging the terms in ( jig) and defining 

a/37,5 K U + ti a ){Uo + h b ) 

ab ^ s (U + E a )(U + E b ) 
x (m* ab m*p - m^m* ) (m Sa m l0 - m so m ia ) , (22b) 

we obtain 



^i 2) (U) = 

Zik Zkj Zjk fki 

^ , , ... f — „ x — " ' n, ,',,/, — ; „ x — ■ '„„/' „'( 

a,b^0 L 



Q ab Tr{a ■ S, £ • S, £ ^™ } 



Zik Z k -i-i k fki 



(23) 



Note that the sum over state labels in Eqs. Q22a| ) and ( |22b| ) is restricted by the conditions 
of Eqs. ( |20a| ) and (|20 b| ) . One notes that the first term in Eq. (^) has a structure similar to 
that of H.^(i,j) in Eq. (|15D and therefore has the same magnetic symmetries. The second 
term, however, leads to magnetic anisotropy even for tetragonal symmetry, for which the 
effective hopping between the coppers is spin-independent. This is elaborated upon in the 
next section. 

Finally we consider the processes in channel b, in which the two holes are on the oxygen 
in an intermediate state and thus experience the Coulomb interaction on the oxygen. These 
processes yield 
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•nn'r,,^' tn tn ' tn l tn \ 




{Up + e n + e n ,){U P + 

fill T Cn' ) 



— Tr j<7 • Sifoni^n'O " ' ^J^On'^no} • 



(24) 



To the order we work, the total magnetic interaction between a pair of copper ions is 
given by 



In the next section we examine 7i(i,j) in the tetragonal phase. In the orthorhombic phase 
the situation is more complicated. In principle, one should start by replacing the tetragonal 
crystal field states by those which are calculated in the presence of the orthorhombic dis- 
tortion. A reasonable approximation is to take account of this distortion by considering the 
crystal field states of the octahedron of oxygen ions assuming this octahedron to be rotated 
rigidly away from its orientation in the tetragonal phase. Then, to the extent that the crys- 
tal field is wholly determined by the octahedron of oxygen neighbors, it will be tetragonal 
in the rotated coordinate system fixed by the shell of oxygen neighbors. Then the result 
of Eq. (p5l) can be used. The major complication is that the hopping matrix elements are 
those between rotated orbitals, as will be detailed elsewhere.il 

Some comments concerning the applicability of our results should be made. As discussed 
in YHAE, we believe that the use of perturbation theory is justified for the calculation of the 
anisotropy of the superexchange Hamiltonian. In contrast, for the isotropic terms it has been 
shownil that perturbation theory is not reliable, mainly because the energies of the excited 
states of the oxygen levels are not very large in comparison to t. However, since the spin- 
orbit interaction takes place on the copper ions, this effect is much reduced in the anisotropic 
terms. In addition, there are contributions to the isotropic interactions which can not be 
obtained by consideration of a single bond. In one of these, a hole hops from a copper to 



(25) 
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one nearest-neighboring oxygen, then diagonally to a different nearest-neighbor oxygen, and 
finally back to the original copper.ll Also, the effective spin Hamiltonian acquires isotropic 
contributions from processes of order t 8 which involve a hole hopping around a plaquette of 
Cu ions.0 These plaquette interactions involve both two-spin interactions (between nearest 
and next-nearest neighbors) and four-spin interactions. 



IV. TETRAGONAL SYMMETRY 

Here we study the effective magnetic Hamiltonian Tl(i,j) [Eqs. (|T5|) and (p3| )-(p5l)l for 
the case of tetragonal symmetry. For this symmetry the effective interaction which is bilinear 
in the spin operators must be of the form 

H(i,j) = J2-J^j)S»(t)S,(j) , (26) 

where \i — x,y, z. Since we are interested in the anisotropic part of the exchange interaction, 
we will drop any contributions which we identify as being isotropic. Obviously, our results 
can not be used for the magnitude of an individual J w , but rather apply to the difference 
between two such quantities. 

Investigation of Eqs. (|15), (p3|) and (p4|) reveals that one can define an effective hopping 



between the copper ions, which is generally given by the product t^t 1 ^^. (The interaction 
Tib(i,j) [Eq. (|24]) 1 requires also the cases i = j.) We therefore start by examining this 



quantity. To this end we note that the nonzero hopping matrix elements between the 
tetragonal states | a) on the coppers and the states | n) on the oxygen are to Px , t\ Px , t yPz 
and t zpy for a bond along the x-direction in the Cu0 2 plane, and analogously t 0py , ti Py , t xpz 
and t zpx for a bond along y.il Using now Eq. (|6]) we find 

taJn'b = tantn'pm* aa mp b 0-(a)0-((3), (27) 

a/3 

where we have omitted the site indices for convenience. It therefore follows that in the case 
n = n' the quantity in Eq. Q2"7p is proportional to the 2x2 unit matrix, namely, the effective 
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hopping between the copper ions is not accompanied by spin-flip. This means that in the 
expressions for the interactions H {l) {i,j) and H^{i,j) [Eqs. (fSJ) and (EJ)] we may take all 
fs outside the trace. Consequently, the contribution to the magnetic energy from Eq. ( |T5D 
is isotropic, and so is that from the first term in Eq. (p3|), proportional to Q a b- The second 



term in Tt^(i,j), which arises from s , leads to an anisotropic magnetic interaction of 
the form of Eq. (|26|) , with 

^=4E^E^Et^- ( 28 ) 

We now consider the magnetic symmetry of Ti,^ b 2 \i,j), which results from the on-site 
Coulomb potential on the oxygen. The only nonzero matrix elements AU nin ' nn i are@ 

AU nn i n i n , AU nn i nn i, AU nnn i n i. (29) 

The first of these implies terms of the form to n i n o [cf. Eq. Q23|)], which are proportional to 
the unit matrix. As a result the first term in the square brackets of (El) disappears. The 
second, which yields an isotropic interaction, can be combined into H.^(i,j) by redefining 
g nn , [Eq. (0)] to be 

2 1 

9nn' Yt 

Uo £n e n' 

1 1 1 2 

Up + AU nn ' n i n + e n + 

^n' ^ ^n ( n' 

For simplicity we have put the AU in the denominator of this expression. The results of 
this paper are correct to first order in A7i c . 

The other two matrix elements of Eq. (EJ]) lead to 

/ 1 1x2/ 1 x2 

AU„ 



/ 1_ J_\2( 1 



E 

nn 1 

x (Tr{<? • Sit rA'o}Tr{<7 • S/on'^no} 

— Trjcr ■ Sit 0n t n r a ■ Sjt 0n /t n0 ^ 
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(Up + 2e n )(U P + 2e n / 

x (^{^ ' Sjton^n'o}Tr|c? • Sjfon^n'o} 
— Tr|(J ■ Sitontn'oV ■ Sjton^n'o}) 



(31) 



The main point to note here is that although ton^n'o implies a spin-flip in the sense that 
this product is proportional to one of the Pauli matrices [cf. Eq. (|27|) 1, it is the same Pauli 
matrix for both products appearing in each term of Eq. (|3l|). This makes the resulting 
interaction isotropic. For example, writing t 0n t n i oc cr M , with \i — x, y or z, we find that the 
spin dependence in each of the terms in Eq. (|31~1) is proportional to 

Tr j<7 • SiCr^jTr j<j • S^} 
- Tr|a • SiU^a • Sjcr^} = 2Sj • S 3 -, (32) 

where we have used the identity (0]). Thus the perturbation processes for which there 
are two holes on the oxygen in the intermediate state do not contribute to the magnetic 
anisotropy in tetragonal symmetry. This result holds because the spin-orbit interaction on 
the oxygen has been discarded. 

The anisotropic exchange interactions in the tetragonal phase thus come exclusively from 
channel a. Also, as previously noted, only the term proportional to K in Eq. (|23|) gives rise 
to anisotropy. Furthermore one notes that there is no anisotropic contribution to J w from 

(2) 

the matrix elements U^ s for which a = (3 or 7 = 5. 

Up to now, our results have been valid to all orders in the spin-orbit coupling, A. We 
now obtain an explicit expression for the anisotropic part of J MM to leading order in A, which 
turns out to be 0(\ 2 ). To this end we use Eq. (0), from which we get 

{1 for a = a 

+ 0(\ 2 ), (33) 
-^ aa - for a ^ a 

where C aa for a 7^ a [see Eqs. (§)] is of order A. Next we note that the summation indices 
a and b of Eq. (p8[ ) can take the values | a),\ b) =\ 1) or | a), \ b) =\ x),\ y) or | z). We 
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therefore need to evaluate J2 n tont n b/^n for | b) =\ 1) and for | b) —\ u), where v = x,y,z. 
Using Eqs. (O) and (|33|) we find 



E 




n 


£n 


E 




n 


£n 


+ E 

n 


tontnl 



E 



= E 

n 



n c " 
+ 0n 



j- 2 _ j.2 /* 



+ 0(A 2 

2 r 



where the fs are the hopping matrix elements for the tetragonal states. Equations 



(34) 
and 



(3|), in conjunction with Eq. ( [22 b|) , imply that the contributions of s where more than 



two of the indices a, (3, 7, and 5 take the values x, y or z are at least of order A 3 . Take 



for example the element u[^ xy . From Eqs. ( 22b ) and (^) we see that the m-products in 



are of order A 2 provided that | b) = | 1) and | a) =| x) or | y). In all other cases the 
m-products are at least of order A 3 . But with this choice [see Eqs. (B^l)] the hopping matrix 



elements provide another factor of A, to render J„„ to be of order A 3 . It follows that out of 



all nonzero U^g^g, the ones that contribute to Jnn up to order A 2 are 



a/375 



(2) 
OizOii 

(2) 



r/ (2) r/ (2) 



U, 



(2) 



tj^s _ rr(2) 



u, 



(2) 



/zl/zO 



U- 



(2) 



1/zO/i' 



(35) 



We are now in position to calculate K% b . When both | a) and | b) are equal to | 1) one finds 

2 
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(U + e x ] 



(61 



^M)A)M rr (2) 
3 ^l^l/i 



£l/x£/xO + £>0n£>fj,l TT (2) 
7 \ ^0/zl/z 



(36) 



When I a) =| 1) and | b) —\ /j) or vice versa, we have 



ft-/* 



Ml/ 



(C/ + ei)(C/o + e M ) 



x 



ei - e M rr - 



£ fO rr (2) 

^O/ilAt 



(37) 



where we have kept terms up to order A since the hopping elements in this case will contribute 
another factor of A. Finally, the case where both | a) and | b) are equal to | fi) requires 
terms to order A and is therefore 
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zU, 



(2) 



(u + e,y w 



(38) 



Combining the results of Eqs. 



-(|38|) we obtain 



J, 



8U< 



(2) 



OpOp, 



r f 2 — f 2 



E 



8U 



(2) 



tontnl 



£p0 y-v 



(C/o + ei) 2 



ton^ 



ol 



16C7, 



(2) 
Oplp 



C 



(C/o + ei)(C/o + e„) 
^ r f 2 _ j.2 

El 71 n 



E 



Un 



£n 

E 



~H n *-n "-^0 ~t~ ^1 n 6n 

It remains to insert here the explicit expressions for the hopping terms, Eqs 
conjunction with the expressions for the angular momentum matrix elements, 
(§). To be specific, we consider a bond along x, and use the values of U a ^& = ( 
listed in Table II. One then finds 
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An alternative derivation of these results is given in the Appendix. 

One should keep in mind that Eq. ( p5|) holds for a single bond. To obtain the effective 
magnetic Hamiltonian of the entire Cu02 plane, one has to sum the magnetic interaction 
7i(i,j) over all bonds, allowing for the crystal symmetry. Within a classical approximation 
the resulting exchange Hamiltonian of the crystal has only an easy-plane anisotropy. To 
obtain the four-fold anisotropy within the easy plane requires a consideration of the spin- 
wave zero point motion.i 

Now we evaluate these results numerically. In this connection it is useful to emphasize 
that for their less general model YHAE have shown that the perturbative results for the 
anisotropy in the J's agrees to within about 10% with the numerical evaluations of exact 
diagonalization within a Cu-O-Cu cluster. For our numerical evaluation we use the param- 
eters of Table 1. We note that all the nonzero hopping matrix elements can be expressed in 
terms of (pda) and (pdn) « — |(pdcr)@: 

to, Px = -V3t 1)Px = ^-{pda) 

t y , Pz = t 0>Py = (pdn) . (41) 

Thereby we find (in /ieV) 

AJ = J av - J zz = 30 , 5J = J ± ~ J|| = 41 . (42) 

where J av = (J xx + J yy )/2, J± = J yy , and Ji = J xx . Note that the gap in the spin-wave 
spectrum due to the easy-plane anisotropy is proportional to (AJ) 1//2 whereas, as explained 
by YHAE, the gap due to the anisotropy within the basal plane is proportional to 5J. 
(This statement applies to systems like YB2CU3O6. In S^CuC^C^ the in-plane anisotropy 
will have contributions from dipolar interactions.) In YHAE, the terms in J xx and J yy 
involving B\/3 were not included because only Coulomb matrix elements involving at most 
two orbitals were kept. To see the effect of the additional terms in the present work, we 
give, for comparison, the perturbative results of YHAE: AJ = 30/ieV and 5 J = 26/zeV. 
(The results from exact diagonalization on a Cu-O-Cu cluster were AJ = 31/ieV and 
5 J = 31/ieV.) 
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V. CONCLUSIONS 



In view of the results of Ref. [|, YHAE already demonstrated that the out-of-plane 
anisotropy in the superexchange interaction between Cu ions is dominated by Coulomb ex- 
change terms. Therefore, the simplified Coulomb interaction of Eq. (|9]), used widely in the 
literature, is insufficient to explain this anisotropy YHAE then considered the anisotropy 
due to the simplest exchange terms, like U a pp a) which involve only two orbitals, and thereby 
obtained an anisotropy in the superexchange interaction whose value agreed with experi- 
ments. 

In the present paper we included all the Coulomb terms which are allowed by tetragonal 
site symmetry, and found that the additional terms, involving (a0| ad), practically do not 
affect the out-of-plane anisotropy A J. The in-plane single-bond anisotropy, 5 J, is some- 
what larger than before. YHAE showed that the in-plane gap in the spin- wave spectrum 
is a manifestation of quantum zero-point fluctuations and is proportional to 5 J. Using the 
results of our present calculation in Eq. (82) of YHAE we estimate the in-plane gap in the 
spin-wave spectrum to be about 33 /xeV ~ 0.27cm -1 . The direct observation of this gap 
would provide an interesting and significant test of our calculations. 
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Foundation. ABH was also supported in part by the U. S. Israel Education Foundation and 
by the National Science Foundation under Grant No. DMR-91-22784. 

APPENDIX A: ALTERNATIVE CALCULATION OF THE EXCHANGE 

ANISOTROPY 

In this appendix we obtain the result for J MA1 for tetragonal symmetry by direct applica- 
tion of perturbation theory. We take the unperturbed Hamiltonian to be the same as that 
of Eq. (^|), except that we do not include in it the spin-orbit interaction on the copper ions. 
[See Eq. (0).] That term is now to be included into 7i\. For simplicity we only consider the 
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contribution to the anisotropy from channel a. A similar calculation to that given here shows 
that channel b gives no anisotropy. Since we only consider channel a, we can eliminate the 
oxygens completely from the problem by introducing an effective hopping matrix element 
t a p between crystal field states on near-neighboring copper ions. For instance, 

^00 = *0,p x / ' = t0,Pxtl,Px/ € Px 1 

t X x — , tyy = t yPz /e Pz , i zz = t zpy /e Py , (Al) 

as in Eq. (53) of YHAE. Accordingly we need to work to fifth order in perturbation theory, 
where the perturbations are V, the spin-orbit interaction, which we take to second order, 
Hhop = T i the hopping between copper ions, which we take to second order, and ATCc — C 
on the copper ions, which we take to first order. We also use the fact that Uap^s is only 
nonzero when a(a)a(fl) ~ a(j)a(5). 

In fifth order perturbation theory there are, a priori, 30 different ways to order these 
perturbations. But obviously the Coulomb perturbation can only appear when the two 
holes are on the same site. So we have only to consider how to insert two powers of V into 
the sequence TCT. There are basically three types of terms to consider. The first is 

HlS { v i r+T ¥)¥i{ v i T+T ¥) • (A2) 

where £ denotes the appropriate energy denominator. The second type of terms are those 
in which the two V's are both to the left of C, and the third type are terms which are the 
Hermitian conjugates of the second type, i. e. those in which the two V's are to the right 
of C. 

That the second type of term does not lead to any anisotropy can be established by the 
following argument. Suppose we show that these terms vanish when applied to any triplet 
spin state. That would imply that these terms are of the form (1/4) — Sj- Sj, which obviously 
gives rise to no anisotropy. The second type of term has first (reading from right to left) 
T acting on the triplet state. That will put the two holes, which were initially distributed 
one on site i, the other on site j, onto the same site, perforce one in state | 0), the other 
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in state |1). Now apply the Coulomb perturbation. This operator can leave the two holes 
in the same states, viz. one in | 0) and the other in In fact, by our observation on the 
form of U a p-ys and by the fact that parallel spins are not allowed in the same spatial orbital, 
there are no other final states. But such a diagonal matrix element of U was treated by 
YHAE and found to give no anisotropy in the second type of term. So we conclude that the 
second type of term does not give rise to any anisotropy. Terms of the third type are the 
Hermitian conjugate of type two and therefore are subject to the same argument. Thus all 
the anisotropic terms are contained in the expression in Eq. ( |A2|) . 



If Tij denotes hopping from i to j, we can write Eq. 



#1 = 2 E Ql^c±Q a , 

a=x,y,z O 



as 



(A3) 



where 



Qa ( Va £ Tij ~\- Tij ^ V a 



(A4) 



Here V a = J2hT a {h)s a (h), where the sum is over the two holes, h. Acting on the ground 
state, the operator Q Q produces two final states, depending on whether the orbital |0) or 
|1) is occupied. So we define 
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C 



77 



a"r\"\a l rj l 



(0 1 diy a r> d iav » AH c 4 av , 4<y<T> I °) 
8 v ",r,'5*'a" (aj'\An c \aj) 



(A7) 



20 



Now we introduce the notation for direct products 



so that 



C' 7 = \0\{a 1 '\m. c \ai) = [0](a 7 | aV) 

in the notation of Eqs. (|). Also [C] = [J J — a • a]/2. 
From Appendix H of YHAE we also take the results 

Q(0)1 = (<- + c 2 ) [Xa a ] - C 2 [Xa a ] [O] 



and 



Q« = C 3 [Xa Q ] 
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and 
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' ex + e a + U S 
e a {ei + U Q ) 



Thus we obtain the result 
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Now we use the equality 



[la a ] [O] [Ia a ] = - [11 - 2a a a a + a ■ a] . 



(A15) 



So, dropping isotropic terms, we hav< 
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Now use 



x \a a a a ] [C x ll + C 2 11 - C 2 0] 
2{al\al) 



A 2 



2(a0\ al) 
AqqAiq, 



■CI [a a a a ] C 3 

C* 3 [a a a a ] [C x ll + C 2 11 - C 2 0] 



2 i al| A a0) [C* x ll + C;il - C*0] [a a a a ] C 3 



(A16) 



2 [a a a a ] [O] = 2 [O] [a a a a ] = [a ■ a - 11] 



(A17) 



Thus the terms involving O are isotropic. So the anisotropic terms are correctly given by 
Eq. (PpT) with 



J, 



~~ A2 l^l + ^2| 
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Now we use Eqs. ( |A12j ), and ( |A13| ) to write 
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This result reproduces Eqs. fl40|) of the text. 



(A18) 



(A19) 
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TABLES 

TABLE I. Values (in eV) of the parameters used to calculate the anisotropic exchange. For a 
discussion of these values, see YHAE. 



A A B C (pda) e t e x = e y e z e Px e Py e Pz 

0.1 7.0 0.15 0.58 1.5 1.8 1.8 1.8 3.25 3.25 3.25 



TABLE II. Values of (a5\(3j) = Ua^s in terms of the Racah coefficients, takenEil from Table 
A26 of Ref. 

a = x a = y a = z 

(q0|q 0) = 3B+C 3B+C C 

(aO\a 1) = -v^B ^3 B 

(al\a 1) = B+C B+C 4B+C 
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